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Abstract 

We reanalyze deep inelastic scattering data of BCDMS Collaboration by including 
proper cuts of ranges with large systematic errors. We perform also the fits of 
high statistic deep inelastic scattering data of BCDMS, SLAC, NM and BFP Col- 
laborations taking the data separately and in combined way and find good agree- 
ment between these analyses. We extract the values of both the QCD coupling 
constant as(M|) up to NLO level and of the power corrections to the structure 
function F2. The fits of the combined data for the nonsinglet part of the struc- 
ture function F2 predict the coupling constant value a;s(M|) = 0.1174 ± 0.0007 
(stat) ±0.0019 (syst) ±0.0010 (normalization) (or QCD parameter = (204 ± 
25 (total exper.err.))MeV^). The fits of the combined data for both: the non- 
singlet part and the singlet one, lead to the values as(M|) = 0.1177 ± 0.0007 
(stat) ±0.0021 (syst) ±0.0009 (normalization) (or QCD parameter = (208 ± 
27 (total exper.err.))Mey). Both above values are in very good agreement with each 
other. We estimate theoretical uncertainties for a<((M|) as ±0.0047 and —0.0057 
from fits of the combine data, when complete singlet and nonsinglet evolution 
is taken into account. 
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1 Introduction 



The deep inelastic scattering (DIS) leptons on hadrons is the basical process to study the 
values of the parton distribution functions which are universal (after choosing of factorization 
and renormalization schemes) and can be used in other processes. The accuracy of the present 
data for deep inelastic structure functions (SF) reached the level at which the Q^-dependence 
of logarithmic QCD-motivated and power-like ones may be studied separately (for a review, 
see the recent papers [|^, and references therein). 

In the present article we analyze at the next-to-leading (NLO) order f\ of perturbative 
QCD the most known DIS SF F2{x,Q^) taking into account SLAG, NMC, BCDMS and BFP 
experimental data We stress the power-like effects, so-called twist-4 (i.e. ~ 



contributions. To our purposes we represent the SF F2{x, Q"^) as the contribution of the leading 
twist part F|'*^'"^(x, Q^) described by perturbative QCD and the nonperturbative part (twist- 
four terms ~ 1/Q'^): 

F,(x, Q') ^ Ft\x^ Q') = F^^^^ix, Q') + ^ (1) 

The SF F2^^^{x, Q^) obeys the (leading twist) perturbative QCD dynamics including the 
target mass corrections (TMC) (and coincides with F|"'^(a;, Q^) when the target mass correc- 
tions are withdrawn). 

The Eq.(|I|) allows us to separate pure kinematical power corrections, i.e. TMC, so that the 
function /i4(x) corresponds to "dynamical" contribution of the twist-four operators. The pa- 
rameterization (|l]) implies that the anomalous dimensions of the twist-four operators are equal 
to zero, that is not correct in principle. Moreover, there are estimations of these anomalous 
dimensions (see fll]]). Meanwhile, in view of limited precision of the data, the approximation 
(P) and one in the footnote 2 give rather good predictions (see discussions in O, |13[). 



Contrary to standard fits (see, for example, O, |T5[) when the direct numerical calculations 



based on Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [16] are used to evalu- 
ate structure functions, we use the exact solution of DGLAP equation for the Mellin moments 



X 



' F2ix,Q'^)dx (hereafter k = fuU,pQCD,tw2,...) 



(2) 



and the subsequent reproduction of F2^""(x,Q^) 



F|"3^^(x,Q2) and/or Fi^\x,Q^) at every 



needed Q^-value with help of the Jacobi Polynomial expansion method [|T7|-||TP[ (see similar 
analyses at the NLO level [T^-[^ and at the NNLO level and above 
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17, 


18 



details in Refs.|19]. Here we consider only some basical definitions in Section 3. 



The paper has the following structure: in Section 2 we present basic formulae, which are 
needed in our analyses: we consider different types of Q^-dependence of SF moments, effects 
of nuclear corrections and heavy quark thresholds, the structure of normalization of parton 

-'^The evaluation of al{Q^) corrections to anomalous dimensions of Wilson operators, that will be done in 
nearest future by Vermaseren and his coauthors (see discussions in [^), gives a possibility to apply many 
modern programs to perform fits of data at next-next-to- leading order (NNLO) of perturbative theory (see 
detail discussions in Summary). 

^The r.h.s. of the Eq.(|l|) is represented sometimes as F2^'^^ {x,Q^){l + hi{x)/Q'^). It implies that the 
anomalous dimensions of the twist- two and twist- four operators are equal to each other. 



densities in singlet and nonsinglet channels. In Section 3 we introduce the basic elements of our 
fits. Sections 4 and 5 contain conditions and results of several types of fits with the nonsinglet 
and singlet evolutions for different sets of data. In Section 6 we study the dependence of the 
results on choice of factorization and renormalization scales. In Section 7 we summarize the 
basic observations following from the fits and discuss possible future extensions of the analyses. 



2 dependence of SF and their moments 

In this section we analyze Eq.(|I|) in detail, considering separately different types of Q^-dependence 
of structure function F^- 



2.1 The leading- twist dependence 

To study the Q^.^ependence of the SF F^^'^[x,Q'^) = F^^{x,Q^) + Fi{x,Q^), which splits 
explicitly into the nonsinglet (NS) part and the singlet (S) one, it is very useful to introduce 
parton distribution functions (PDF)[|: gluon one fdx, Q^) and singlet and nonsinglet quarks 
ones fsix,Q'^) and fNsix,Q'^). 

The moments M^^(Q2) and M^{Q^) of NS and S parts of SF F2 (see Eq.(|) for definition) 
are connected with the corresponding moments of PDF fi{x, Q"^) (hereafter i = NS, S, G) 

Mn,Q')= [\^-^fi{n,Q')dx 
Jo 

in the following way (see [^], for example) 

M^'iQ') = K^s{f)-CTs{nMQ'))-fNs{n,Q^) 



(3) 



Ks{f) ■ Cf\n,as{Q')) ■ fs{n,Q') + C'j'{n,UQ')) " fG{n,Q 



whereQ 



(4) 



and Cj'^'^{n,as{Q^)) are so-called Wilson coefficient functions. We have introduced here also 
the coefficients 



/ 



Ksif) = E^-//' ^^sif) = el-Ksif), 



(5) 



m=l 



which come from definition of SF F2 (see, for example, |^). Here / is the number of active 
quarks and is charge square of the active quark of m flavor. 

2.1.1 The dependence of as(Q^) is given by the renormalization group equation, which in 
the NLO QCD approximation reads: 



1 



1 



+ ^ 



In 



MQ') (/9o + Aa.(M|)) 



(6) 



"^Our PDF are multiplied by x to compare with standard definition. 

^Sometimes the last term Ks{,f) ■ CQ'^{n,as{Q^)) ■ foin^Q^) we will call as gluon part of singlet moment 



and denote it as M^{Q'^). 



where as(M|) is served as a normahzation. Here and below we use and f3i for the first and 
the second terms with respect to of QCD /3-function: 



- (3oal - (3{al + 



The equation (|^) allows us to eliminate the QCD parameter Aqcd from our analysis. How- 
ever, sometimes we will present it in our discussions, essentially to compare it with the results 
of old fits. The coupling constant as(Q^) is expressed through Aqcd (in MS scheme, where 



A 



QCD 



A 



MS. 



as 



(/?o + /3ia.(g2)) 



AI- 
MS, 



(7) 



The relation between the normalization as(M|) and the QCD parameter Aqcd can be 
obtained from Eq.(^ with the replacement M|. 

We would like to note that the approximations of Eq. (|^) , based on the expansion of inverse 
powers of In (Q^/A^^) are very popular. The accuracy of these expansions for evolution of Ug 
from O(GeV^) to M| may be as large as 0.001 |]13|, which is comparable with the experimental 
uncertainties of the as(M|) value extracted from the data (see our analyses in Sections 4 and 
5). 

Note also that sometimes (see, for example, E^) the equation 



(8) 



as{Q 

is used in analyses. This equation can be obtained from the basic equation 

Q2 >^ i^MQ') db 



In 



A^ 

MS, 



(9) 



by expansion of inverse QCD /3-function in r.h.s. of (^ l//3(as) in powers of a^. The difference 
between Eqs.(^) and (|^) may be as large as 0.001 at O(GeV^) range. In order to escape the 
above uncertainties we use in the analyses the exact numerical solution (with accuracy about 
10~^) of Eq.(|D instead. For recalculation of the QCD parameter Aj^ from A^ to A^^^ (i.e. 
from A^ at / active quark flavors to Aj^ at / ± 1 active quark flavors), because Po and /5i are 
/-dependent functions, we use formulae at NLO approximation from Ref. ||29[ (see discussions 
in the subsection 2.4). 



2.1.2 The coefficient functions Cf {n, as{Q )) {i = NS, S, G) have the following form 



Cf^(n, a,(Q")) = 1 - 5f + a, ■ Bi{n) + O(a^) (5™ is Kronecker symbol). 



(10) 



where the NLO coefficients Bi{n) are exactly known (see, for example, p^ ). 

The Q^-evolution of the moments fi{n.,Q'^) is given by the well known perturbative QCD 
28| , |30| formulae: 



fNs{n,Q^) 
fNs{n, Ql) 



asm 



7^"^(n)/2/3o 



fA^,Q') = fnn,Q') + fiin,Q') ij = S,G) 



(11) 



ftMl) 



where^ 



7fH 



efi{n)fi{n,Q'), = S,G) 

1 



(12) 



2 



(0) 
SG 



2V -(0)^--^ -(0)^--^ ^ J''^^ 



r± (n) - (n) 



(0), 



(j ^ 



The functions H^^{n, Q"^, Ql) and Hp{n, Q^) 

are nonzero above the leading order (LO) 

approximation and may be represented as 



H'^'in^Q^Ql) = 1+ (a,iQ')-a,iQl)^Z^sin) + 0(aliQ') 
Hf{n,Q\Ql) = l+{as{Q^)-as{Ql))z±±{n) 



(13) 



where 



+ 

ZNs{n) 
Z±±in) 



as{Ql 



a.(Q2). 



2Po 
1 



(1) /X (0) / X A 

^r±±{n)-r±'{n)■ — 



7g( 



n) 



and 



7S(n) 



zU^) 



Z^n) 



'GG 
'55 



7S( 



n] 



j,l=S,G 

^jG^Gi ~ ^55 755 + (^G5 ~ ^GG/^Scj^GG 



(14) 
(15) 

(16) 
(17) 



j=S,G 
„(0) 



As usually, here we use 7]v5(^), 7]°^^) (j, ^ = S,G) and 7]v5('^)' 7^7^'' l'^) ^^e first and 
the second terms with respect to of anomalous dimensions 77V5(^;^s) and 7j;(n, a^) (see, for 
example, p5[]). 



(1) 



2.1.3 In this subsection, we would like to discuss a possible dependence of our results on the 
factorization scale np and the renormalization scale /ir, which appear (see, for example, [p^, |36|1) 
because perturbative series are truncated. These scales fip = kpQ^ and fij^ = kji^p = kjikpQ^ 
can be added to the r.h.s. of the equations (||) and (pA]), respectively. 

Then, the equations (|^) are replaced by 



KNsif) ■ C'j:i's\n,asikpQ')) ■ hsin, kpQ' 



KsU) ■ Cf\n,as{kpQ')) ■ fs{n, kpQ') + C'^\n,a,{kpQ')) ■ fain, kpQ 



^We use a non-standard definition (see [^if ) of the projectors e^(n), which is very convenient beyond LO 
(see Eq. (|l7| ) and Q). The connection with the more usual definition a, a and e in ref. jsj, ^ is given 
by: eaa{n) = a{n), eanin) = a(n) and eQa{n) = e{n) 



The equations (0) are replaced correspondingly by 



fNsjn.kpQ'^) 



asjkFkRQ. 
asikpkRQ^ 

asjkpkRQl 
ttsikEkuQ"^ 



■H^'{n, kpknQ'',kFknQ, 



(19) 



7^ (n)/2/3o 



The coefficients Cms. Cs, Cg, and Hf can be obtained from the ones Cj^s, Cs, Cq, 



H^^ and by modification in the r.h.s. of equations (p!0D, ( p!4D and ([T5|) as follows: 



inEq. (0) 



BNsin) 



asikpQ'^ 



BNs{n) + |7{?sH ^^^F 



B,{n) 



7gHlnfc^ {j = S,G) 



in Eqs. (|T4D and (|T|) 



Z±±(n) 



2 '^-^ 



asikpkRQ^ 



ZNs{n) + ^'y^^l{n)\nkii 



(20) 



(21) 



(22) 



(23) 



The Eqs. (|2l| ) can be obtained easily using, for example, the results of The Eqs. (^) 
can be found from the expansion of the coupling constant asikpkRQ"^) around the one cisiQ^) 
in the r.h.s. of the exact solution of DGLAP equations (see Eqs.(|19]) and (p!3|)). 

The changes ( pi]) and ( p3|) of the results for Q^-dependence under variation of fc^ and 
kpi (usually from 1/2 to 2) give an estimation of the errors due to renormalization and 
factorization scale uncertainties. Evidently that, by definition, these uncertainties are connected 
with the impact of unaccounted terms of the perturbative series and can represent theoretical 
uncertainties in values of ffited variables. Indeed, an incorporation of NNLO corrections to the 



analysis strongly suppress these uncertainties (see ||39|, |iO|). 

We study exactly the and iir dependences here for fitted values of couphng constant. 
The results of the study are given in the Section 6. 

As one can see in Eqs. (^) and (p2|), the coupling constant has different arguments in the 
NLO corrections of coefficient functions Cns and Cj (j = S", G) and in the NLO corrections H^^ 
and H"^ of the Q^-evolution of parton distributions. We would like to note that the difference 
between the corresponding coupling constants TisikpQ'^) and as^kpkRQ"^) is proportional to 
and, thus, mathematically negligible in our NLO approximation. 

Then, we can use the replacement ( ^2|) in coefficient functions too, as it has been done in 



previous studies [|38|, ^ We note that the replacement asi^kpQ'^) — >■ agikpkRQ'^) in Eq. 

(pO|) increases slightly the factorization scheme dependence of the results for coupling constant 
(see analyses based on nonsinglet evolution and discussions in Section 6). 



^In the recent articles ^ ^ the variation from 1/4 to 4 has been used. In our opinion, the case 
kp = kii — A leads to very small scale of coupling constant: (5^/16, that requires to reject many of experimental 
points of used data, because we have the general cut > 1 GcV^. So, we prefer to use the variation of scales 
from 1/2 to 2. 



2.2 Normalization of parton distributions 

The moments fi{n, Q'^) at some Qq is theoretical input of our analysis which is fixed as follows. 

For fits of data at x > 0.25 we can work only with the nonsinglet parton density and use 
directly its normalization fNs{x-,Ql) (see, for example, P^-[PB|]): 

fNs{n,Ql) = / dxx'^'^fNs{x,Ql), 
Jo 

fNs{x,Ql) = A^s{Ql)il-x)''-'^^?^\l + dNsiQl)x) (24) 
where Ans{Qo), ^NsiQo) dNsiQl) some coefficients f\. 

At the analyses at arbitrary values of x we should introduce the normalizations for densities 
of individual quarks (g = u,d,s,...) and antiquarks (g = u,d,s,...) fg{x,Ql) and fg{x,Ql) 
having the moments: 

Mn,Ql) = ['dxx--'f,{x,Ql), (25) 

The distributions of u and d quarks fu{x,Ql) = u{x,Ql) and fd{x,Ql) = d{x,Ql) are 
split in two components: the valent one Uv{x,Ql) and dy{x,Ql) and the sea one Uv{x,Ql) 
and dy{x,Ql). For other quark distributions and antiquark densities we keep only sea parts. 
Moreover, following [^, we suppose equality of all sea parts and mark their sum as S{x, Qq). 

We use the following parameterizations for densities m„(x, Ql), Uy{x, Ql), S{x, Ql) fcix, Qq): 

B{a,„(Qi), 0„(Q5) + 1) 



f] (r O'^) = = jMQo) (I - rf-'^^o) (26) 

Mx,Qo) Bia,{Ql)MQl) + l) ^ ^ ' ^ ' 

S{x,Ql) = Cs{Ql)x^'^'^o)^i-:,fs(Ql)^ 

fGix,Ql) = Cc(Q^)x'^«(«o)(i_^)^o(Qg)^ (27) 



where B{a, b) is the Euler beta-function. The parameterizations ( p6D have been chosen to 
satisfy (at the normalization point Qq) the known rule: 



1 

dxV{x,Q^) = 3, 



where V{x, Q^) = Uy{x, Q^) + dy{x, Q^) is the distribution of valent quarks. 

We note that the nonsinglet and singlet parts of quark distributions fNsix, Qq) and fsix, Qq) 
can be represented as combination of quark ones 

/ 

fs{x,Ql) ^ 5^/,(x,Q^) = V{x,Ql)+Six,Ql) (28) 
fNs{x,Ql) = Uy{x,Ql) - dy{x,Ql), (29) 



^We do not consider here the term ~ ^''"^''^o' in the normalization fNs{x, Qq), because x > 0.25. The correct 
small-x asymptotics of nonsinglet distributions will be obtained by Eq. (|2^) from the corresponding parameters 
of the valent quark distributions (|2^ ) fitted with complete singlet and nonsinglet evolution in Section 5. 



where the r.h.s. of Eq. (p9D is correct only in the framework of our supposition about equahty 
of antiquarks distributions and sea components of quark ones. 

In principle, following the PDF models used in [|l^, |12| and above Eq.(^) one can add in 
Eq.(p7D terms proportional to \fx and x. However, the terms ~ ^Jx are important only in the 
region of rather small x (see discussion in The terms ~ x lead only to replacement of C^, 

Oj and hi values (see, for example, |^). Thus, we neglect these terms in our analysis. 

In the most our fits we do not take into account also the terms ~ x°°*^'^o^ and ~ x°'^^'^^^ 
into gluon and sea quark distributions, because we do not consider experimental data at small 
values of Bjorken variable x 0. We hope to include HI and ZEUS data [^, Q in our future 
investigations and then to study Q^-dependence of the coefficients aG(Q^) and as((5^). 



which could be very nontrivial (see, for example, Refs. H3, H^, |32| and references therein) 



We impose also the condition for full momentum conservation in the form: 

1 = Pam^rPM). 



(30) 



where 



dxfG{x,Q'^ 



dx[fNs{x,Q^) + fs{x,Q^ 



(31) 



The coefficients Ci{Ql), ai{Ql), bi{Ql), Ci{Ql) and di{Ql) should be found together with 
/i4(x) (see subsection 2.6) and the normalization asi^M^) of QCD coupling constant (or QCD 
parameter A) by the fits of experimental data. 



2.3 Target mass corrections 



The target mass corrections [^, ^ modify the SF F2 {x, Q ) in the following way 



1^ 

7372 2 



<Y>6 



6 



Mlucl. 



+ 12 



MLci. 



X 



d^ T7tw2(c' n2^ 



(32) 



a/1 + x'^M'^^^i IQ^ and the Nachtmann variable 



where Mnud is the mass of the nucleon, r 
i = 2x/(l + r). 

In our analyses below, we will use all this (|3^) representation |. We would like to keep the full 
value of kinematic power corrections, given by nonzero nucleon mass. Then, the excess of 1/Q^ 
dependence encoded in experimental data will give the magnitude of twist-four corrections, 
which is most important part of dynamical power corrections. 



^However, we have performed several fits with nonzero aQ and 05 values taken into account (see Section 
5). We have found a negative value of them: aQ — as ^ —0.18 (that is in agreement with [^) but these 
results cannot be considered seriously without taking into account HI and ZEUS data Q (see, however, 
discussions in the subsection 5.3.4). 

^It is contrary to jl3j, where only the term ^ M^^^i /Q^ has been used. We note that the appearance of the 
terms ~ M^^^^JQ^ at x ^ 



1 (see, for example, |50 ), i.e. the absence of the equality F2*^'^^(l, 



= 0, is not 

important in the our analyses because we do not use experimental data at very large x values: a: < 1. 



2.4 Thresholds of heavy quarks 



Modern estimates performed in [|51| , |52| have revealed a quite significant role of threshold effects 
in the as{Q'^) evolution when the DIS data lie close to threshold points = ~ "^z+i (^^ 

the position of so-called "Euclidean-reflected" threshold of heavy particles). The corresponding 
corrections to the normalization as(M|) can reach several percent, i.e. , they are of the order 
of other uncertainties which should be under control at our analysis. 

An appropriate procedure for the inclusion of threshold effects into the Q^-dependence of 
«s(Q^) in the framework of the massless MS scheme was proposed more than 10 years ago 
| pB| , |M[ : transition from the region with a given number of flavors / described by massless 
asiQ"^', f) f^to the next one with / + 1 ("transition across the Mj+i threshold") is realized here 
with the use of the so-called "matching relation" for as{Q^) 0|. The latter may be considered 



as the continuity condition for as{Q ) on (every) heavy quark mass rrif 



a^iQ' = M^f^,-J) = as{Q^ = M';f + l) and (33) 



Mj+i = m;+i (34) 

that provides an accurate as((5^)-evolution description for Q'^ values not close to the threshold 
region (see and references therein). 

At the analyses based on nonsinglet evolution, the additional /-dependence comes only 
from the NLO correction of NS anomalous dimensions (see [^) |^. In the Section 3 we check 



numerically the dependence of the results from the matching point. We use two matching 
points: (Q) one and 

My+i = 2mj+i (35) 

and demonstrate very little variations of the results Q (see Section 4 and discussion there). 

As we know, for singlet part of evolution no simple recipe exists for exact value of the 
matching point Mf^i. From one side, as in the nonsinglet case, there is Q^-evolution of the SF 
moments which leads to above condition (0). But here we have also the generation of heavy 
quarks (at lowest nontrivial order, in the framework of the photon-gluon fusion process), that 
gives contributions to gluon part of the singlet coefficient function. The photon-gluon fusion 
needs the matching point at the value of Q^, when = 4:m'j_^^, i.e. 

^J+i ■ ^ + Mi,, = 4m5+i (36) 

At small X values the condition (|36| ) is quite close to the one ( pi]) (for example, at x = 0.2 
= — Mnuci.l^)i but at the range of large and intermediate values of x, the 

value of Mj+i is essentially large to compare with one of Eq.(p^. For example, at x = 0.5 
M'jj^^ = 4:m'j_^-^ — M^^^i , that is very close to the matching point (^5]). At larger x values the 
value of Mj^^ will be close to ones in 



We note, that the difference between nonsinglet and singlet Q^-dependences comes from 
contribution of gluon distribution. The contribution is negligible at x > 0.3 that supports 
qualitatively the choice (|3^ ) as the matching point. 



^"Following 1^ in this subsection we use the form as{Q^;f) for the coupUng constant with purpose to 
demonstrate its /-dependence through the ones of /3o and /3i coefficients. 

^^The corresponding moments at any value are proportional to same coefficient K^sif)- Thus, the 
coefficient can be always taken up by the normalization M^^{Q'^). 

^^We will not take into account a small variation (see |Q) of the continuity condition ( p3| ) because of the 
matching point (Ba). 



We would like to note also, that at NLO approximation and above, the situation is even 
more difficult in singlet case, because every subprocess generates itself matching point Mj+i to 
coefficient functions. To estimate a possible effect of a dependence on matching point, we will 
fit data (in Section 3) with two different matching points: (0) one and (|35| ) one. Surprisingly, 
at the singlet case, where all functions coming to Q^-evolution are /-dependent, we do not find 
a strong /-dependence of our results (see Section 5 and discussions there). 



2.5 Nuclear effects 



Starting with EMC discovery in , it is well known about the difference between PDF in free 
hadrons and ones in hadrons in nuclei. We incorporate the difference in our analyses. 

In the nonsinglet case we parameterize the initial PDF in the form ( p4D for every type of 
target. We have 



fNsin,Q't>) = l' dxx--'fU^,Ql), (37) 
where 

fU^,Ql) = KsiQDi^ - + dUQ» (38) 

and A = H, D, C and F in the case of if2, -D2, and Fe^^ targets, respectively. 

In the singlet case we have many parameters in our fits, which should be fitted very carefully. 
The representations similar to ( p8|) for gluon and sea quark PDF should complicate our analyses. 
To overcome the problem, we apply the Eqs. ( pTf ) and (^) only to H2 and D2 cases. For heavier 
targets we apply simpler representations for F2 structure functions in the form: 

F^{x,Q^) = F,^{x,Q')-Kt{l-K^x + K^x') {A = , Fe''), (39) 



where we use experimental observation (see |^0[ and references therein) about approximate 
Q^-independence of EMC ration F^/F^. 

2.6 Higher- twist corrections 

For n-space eq.(|I|) transforms to 

Mf" (Q^) = Mf^^iQ^) + (40) 
where h^ln, Q"^) are the moments of the function h^{x, Q^): 



hiin) = / ^hi{x)dx (41) 





The shape h/^{x) (or coefficients hiin)) of the twist-four corrections are of primary consid- 
eration in our analysis. They can be chosen in the several different forms: 

• the twist-four terms (and twist-six ones) are fixed in agreement with the infrared renor- 
malon (IRR) model (see |I], ^ and references therein). 



^•^The small dependence of EMC ration has been observed also in theoretical studies. For example, in the 
framework of rescaling model the dependence is very small (see |^). It has double-logarithmic form 
and locates only in argument of Euler NP-function. 



The twist-four term in the form h^^x) ~ -^lnF,^^{x,Q'^) ~ 1/(1 — a;) (see and 



references therein). This behavior matches the fact that higher twist effects are usually 
important only at higher x. The twist-four coefficient function has the form Cf'^'^{n) = 
{n - l)Af\ 

• The twist-four term h^lx) is considered as a set free parameters at each bin. The set 
has the form h{^^^{x) = ^[^j^ h^i^i), where / is the number of bins. The constants /i4(xj) 
(one per x-bin) parameterize x-dependence of h^^^^[x). 

The first two cases have been considered already in ||2l| and will be studied carefully later 
45|. Here we will follow the last possibility Q. 



3 Fits of F2'. procedure 

To clear up the importance of HT terms we fit SLAC, NMC, BCDMS and BFP experimental 
data 0-|]lO| (including the it systematic errors), keeping identical form of perturbative part at 
NLO approximation. In the Section we demonstrate the basic ingredients of the analyses. 

As it has been already discussed in the Introduction we use the exact solution of DGLAP 
equation for the Mellin moments M^'^{Q'^) (^ of SF F2^'^{x,Q'^) and the subsequent repro- 
duction of F/""(x,Q2)^ F|^^^(x,Q2) and/or F^'^^x.Q^) at every needed Q^.yalue with help 
of the Jacobi Polynomial expansion method. The method of the Jacobi polynomial expansion 



was developed in |T8| and described in details in Refs. ||T8[| . Here we consider only some 
basical definitions. 

Having the QCD expressions for the Mellin moments M^{Q^) we can reconstruct the SF 
F|(a;,Q2) as 

^771 ax ^ 

F^^''—{x,Q')=x%l-x)'' e:'''(a;)5^cfHa,/3)Mj^2(Q'), (42) 

n=0 j=0 

where 0^' are the Jacobi polynomials Q and a, b are the parameters, fitted by the condition of 
the requirement of the minimization of the error of the reconstruction of the structure functions 
□ (see Ref.pl for details). 



First of all, we choose the cut > I GeV^ in all our studies. For Q"^ < 1 GeV^, the 
applicability of twist expansion is very questionable. 

Secondly, we choose quite large values of the normalization point Q^. There are several 
reasons of this choice: 

• Our above perturbative formulae should be applicable at the value of Q^. Moreover, 
the higher order corrections ~ a'^{Ql) (n > 2), coming from normalization conditions of 
PDF, are less important at higher Ql values. 

• It is necessary to cross heavy quark thresholds less number of time to reach = M|, 
the point of QCD coupling constant normalization. 



^■^In conclusion we present, however, several comments about an application of higher-twist corrections in the 
form of IRR model. 

-'^^We would like to note here that there is similar method based on Bernstein polynomials. The method 
has been used in the analyses at the NLO level in joj, ^ and at the NNLO level in j6^, . 

There is another possibility to fit data. It is possible to transfer experimental information about structure 
functions to their moments and to analyze directly these moments. This approach was very popular in the past 
(see, for example, Ref. |Q) but it is used very rarely at present (see, however, and references therein) 
because a transformation of experimental information to the SF moments is quite a difficult procedure. 



• It is better to have the value of Qq around the middle point of logarithmical range of 
considered values. Then at the case the higher order corrections ~ (as(Q^) — as{Ql))"' 
(n > 2) are less important. 

Basic characteristics of the quality of the fits are x^/DOF for SF F2 and for its slope 
dlnF2 / dlnQ'^ , which has very sensitive perturbative properties (see p8| ). 

As these fits involve many free parameters independent of perturbative QCD, it is important 
to check whether, in the results of the fits, the features most specific to perturbative QCD are 
in good agreement with the data. The slope dlnF2/ dlnQ"^ has really very sensitive perturbative 
properties and will be used (see the Figs. 4-8 and 10-12) to check properties of fits. Indeed, 
the DGLAP equations predict that the logarithmical derivations of SF and PDF logarithms 
are proportional very nearly to coupling constant as(Q^) with an x-dependent proportionality 
coefficient that depends (at x > 0.2) only weakly on the x-dependence of the SF and PDF. 
Thus, the study of the Q^-dependence of the slope dlnF2/dlnQ'^ leads to obtain the direct 
information about the corresponding Q^-dependence of QCD coupling constant and to verify 
the range of accuracy for formulae of perturbative QCD. 



We use MINUIT program |68| for minimization of two x values: 



X\F2 



rpexp 
^2 



rpteor 
^2 



and (slope) 



j~^exp j~^teor 



/^JJexp 



D 



dlnF2 
dlnQ^ 



We would like to apply the following procedure: we study the dependence of x^/DOF value 
on value of cuts for various sets of experimental data. The study will be done for the both 
cases: including higher twists corrections (HTC) and without them. 

We use free normalizations of data for different experiments. For the reference, we use the 
most stable deuterium BCDMS data at the value of energy Eo = 200 GeVQ. Using other types 
of data as reference gives negligible changes in our results. The usage of fixed normalization 
for all data leads to fits with a bit worser x^. 



4 Results of fits of F2: the nonsinglet evolution part 

Firstly, we will consider the Q^-evolution of the SF F2 at the nonsinglet case where there are the 
contributions of quark densities only and, thus, the corresponding fits are essentially simpler. 
The consideration of the nonsinglet part limits the range of data by the cut x > 0.25. At 
smaller x-values the contributions of gluon distribution is not already negligible. 

Hereafter at nonsinglet case of evolution we choose Qq = 90 GeV^ for the BCDMS data and 
the combined all data and Ql = 20 GeV^ for the combined SLAC, NMC, BFP one, respectively. 
The choice of Qq- values is in good agreement with above conditions (see the previous Section). 
We use also Nmax = 8, the cut 0.25 < x < 0.8. The A^maa-dependence of the results has been 
studied carefully in Ref. |TB| (see also below the Table 3). 

4.1 BCDMS C^^ + H2 + D2 data 

We start our analysis with the most precise experimental data [0, |^, ^ obtained by BCDMS 
muon scattering experiment at the high values. The full set of data is 607 points (when 
X > 0.25). The starting point of QCD evolution is Ql = 90 GeV^. 

It is well known that the original analyses given by BCDMS Collaboration itself (see also 
Ref. [0) lead to quite small values of q;s(M|): for example, as(M|) = 0.113 has been obtained 

Eq is the initial energy lepton beam. 



in [T^ Although in some recent papers (see, for example, [|T^, 13, 69|) more higher values 
of as(M|) have been observed, we think that an additional reanalysis of BCDMS data should 
be very useful. 

fsee also [I7l|, 



Based on study 



d9|) we propose that the reason for small values of as{M'^) 
coming from BCDMS data is the existence of the subset of the data having large systematic er- 
rors. Indeed, the original analyses of H2, D2 and C^^ data performed by BCDMS Collaboration 
lead to the following value of QCD mass parameter (see Refs. 0, ||, 



A. 



(4^ 
MS 



220 ± 13 (stat) ± 50 (syst) ) MeV 



(43) 



i.e. the systematic error is four times bigger than the statistical one. Hereafter the symbols 
"stat" and "syst" mark the statistical error and systematic one, respectively. 

We study this subject by introducing several so-called F-cuts (see and subsections 
4.1.1 and 5.1.1). Excluding this set of data with large systematic errors leads to essentially 
larger values of as(M|) and very slow dependence of the values on the concrete choice of the 
F-cut (see below). 



4.1.1. The study of systematics. 



The correlated systematic errors of the data have been studied in |72|, together with the 
other parameters. Regions of data have been identified in which the fits cause large systematic 
shifts of the data points. We would like to exclude these regions from our analyses. 

We have studied influence of the experimental systematic errors on the results of the QCD 
analysis as a function of Ycua, Ycuti and Ycutb applied to the data. We use the following x- 
dependent y-cuts: 



y 


> 0.14 


when 


0.3 < 


X 


< 0.4 


y 


> 0.16 


when 


0.4 < 


X 


< 0.5 


y 


^ Ycut3 


when 


0.5 < 


X 


< 0.6 


y 


> YcuU 


when 


0.6 < 


X 


< 0.7 


y 


> Ycut5 


when 


0.7 < 


X 


< 0.8 



(44) 

We use several sets N of the values for the cuts at 0.5 < x < 0.8, which are given in the 
Table 1. 



Table 1. The values of Y^uts, Y^utA and Fc^ts- 








1 


2 


3 


4 


5 


6 


Ycut3 





0.14 


0.16 


0.16 


0.18 


0.22 


0.23 


YcutA 





0.16 


0.18 


0.20 


0.20 


0.23 


0.24 


y 

-' cuts 





0.20 


0.20 


0.22 


0.22 


0.24 


0.25 



The systematic errors for BCDMS data are given |]^, ^ as multiplicative factors to be 
applied to -^2(0;, Q^): /r, /?,, /«, fd and fh are the uncertainties due to spectrometer resolution. 

We would like to note that the paper has a quite strange result. Authors of the article have obtained 
the value — 263 MeV, that should lead to the value of coupling constant as(M|) is equal to 0.1157. 

^^Hereafter we use the kinematical variable Y = (i?o — E)/Eq, where Eq and E are initial and scattering 
energies of lepton, respectively. 




Figure 1: The study of systimatics at different i/cut values in the fits based on nonsinglet 
evolution. The QCD analysis of BCDMS C^^ , H2, D2 data (nonsinglet case): Xcut = 

0.25, 

Ql = 90 GeV^. Thresholds of c and b quarks are chosen at = 9 GeV^ and = 80 GeV^, 
respectively. The inner (outer) error-bars show statistical (systematic) errors. 



beam momentum, calibration, spectrometer magnetic field calibration, detector inefficiencies 
and energy normalization, respectively. 

For this study each experimental point of the undistorted set was multiplied by a factor 
characterizing a given type of uncertainties and a new (distorted) data set was fitted again in 
agreement with our procedure considered in the previous section. The factors (/r, ft, fs, fd, fh) 
were taken from papers 0, ||, ^ (see CERN preprint versions in ||^, ||, ^j). The absolute differ- 
ences between the values of for the distorted and undistorted sets of data are given in Table 
2 and the Fig. 1 as the total systematic error of as estimated in quadratures. The number of 
the experimental points and the value of ag for the undistorted set of F2 are also given in the 
Table 2 and the Fig. 1. 

Table 2. The values of as(M|) at different values of A^. 



N 


number 


X\F2)/D0F 


as{90 GeV^) ± stat 


full 




of points 






syst. error 





607 


1.03 


0.1590 ± 0.0020 


0.0090 


1 


511 


0.97 


0.1711 ± 0.0027 


0.0075 


2 


502 


0.97 


0.1720 ± 0.0027 


0.0071 


3 


495 


0.97 


0.1723 ± 0.0027 


0.0063 


4 


489 


0.94 


0.1741 ± 0.0027 


0.0061 


5 


458 


0.95 


0.1730 ± 0.0028 


0.0052 


6 


452 


0.95 


0.1737 ± 0.0029 


0.0050 



From the Table 2 and the Fig. 1 we can see that the values are obtained for = 1 ^ 6 
of Ycut3, Ycuti and Ycutb are very stable and statistically consistent. The case = 6 reduces the 



systematic error in by factor 1.8 and increases the value of a^, while increasing the statistical 
error on the 30%. 

After the cuts have been implemented (in this Section below we use the set = 6), we have 
452 points in the analysis. Fitting them in agreement with the same procedure considered in 
the previous Section, we obtain the following results: 

a, (90 GeV^) = 0.1737 ± 0.0029 (stat) ± 0.0050 (syst) ± 0.0025 (norm) 

(45) 

a,(M|) = 0.1153 ± 0.0013 (stat) ± 0.0022 (syst) ± 0.0012 (norm), 

where hereafter the symbol "norm" marks the error of normalization of experimental data. 
Thus, the last error (±0.0011 to as(M|)) comes from difference in fits with free and fixed 
normalizations of BCDMS data 0, H, 3 having different values of energy. 

So, for the fits with NS evolution of BCDMS data |^, ^ with minimization of systematic 
errors, we have the following results: 

a^(M|) = 0.1153 ±0.0028 (total experimental error) (46) 

Here total experimental error is squared root of sum of squares of statistical error, systematic 
one and error of normalization. 

The value of a^j(M|) corresponds to the following value of QCD mass parameter: 

A^^ = (l81± 32 (total experimental error) ) MeV, 

A^ = (257 ± 40 (total experimental error)) MeV (47) 



4.1.2. The study of A^max-dependence. 



Following to |T^, we study the dependence of our results on the Nmax value. The full 
set of data is 452 points. The Q^-evolution starts at Qq=90 GeV^. 

As it can be seen in the Table 3, our results are very stable, that is in very good agreement 



with p 



Table 3. The values of as(M|) at different values of N„ 



^max 


X\F2)/D0F 


(slope) 
for 6 points 


a, (90 GeV^) ± stat 
stat = 0.0038 


as(M|) ± stat 
stat = 0.0013 


3 


1.08 


7.3 


0.1720 


0.1155 


4 


0.97 


11.3 


0.1715 


0.1143 


5 


1.11 


6.9 


0.1729 


0.1144 


6 


0.95 


3.6 


0.1747 


0.1157 


7 


0.94 


5.4 


0.1740 


0.1154 


8 


0.94 


6.8 


0.1738 


0.1153 


9 


0.94 


7.6 


0.1735 


0.1152 


10 


1.07 


7.7 


0.1735 


0.1152 


11 


1.08 


7.2 


0.1726 


0.1149 


12 


1.04 


7.1 


0.1731 


0.1152 


13 


1.11 


7.1 


0.1725 


0.1149 



Starting with N^ax = 5, where our results are aheady very stable, we put the results 
together and can calculate average value of as(M|) = 0.1152 and estimate average deflection. 
The deflection is 0.0002 and can be considered as error of the Jacobi Polynomial expansion 
method, i.e. method error. 

4.2 SLAC and NMC H2 + D2 data and BFP Fe data 

We continue our NS evolution analyses by fits of experimental data @, H, H, [lOl obtained by 
SLAC, NM and BFP Collaborations. The full set of data is 345 points (when x > 0.25): 238 
ones of SLAC, 66 ones of NMC and 41 ones of BFP. The starting point of QCD evolution is 
Ql = 20 GeV^, the Q^.cut is > 1 QeV^. 

For illustration of importance of 1/Q^ corrections we fit the data in the following way. First 
of all, we analyze the data applying only perturbative QCD part of SF F2, i.e. Fg'"^. Later, we 
have added l/Q"^ corrections: firstly, target mass ones and later twist-four ones. As it is possible 
to see in the Table 4, we have the very bad fit, when we work only with twist-two part Fg""^. 
The agreement with the data is improved essentially when target mass corrections have been 
added. The incorporation of twist-four corrections leads to very good fit of the data. Neglect 
of systematic errors deteriorates twice our results. We combine the statistical and systematic 
errors in quadrature. 



Table 4. The values of q;s(M|) and different regimes of fits. 



N of 


TMC 


HTC 


syst. 


X\F2)/D0F 


(slope) 


a, (20 GeV^) 


«s(M|) 


fits 






error 




for 6 points 


± stat 




1 


No 


No 


Yes 


6.0 


1050 


0.2131 ± 0.0012 


0.1167 


2 


Yes 


No 


Yes 


2.3 


224 


0.2017 ± 0.0013 


0.1133 


3 


Yes 


Yes 


No 


1.8 


12.0 


0.2230 ± 0.0030 


0.1195 


4 


Yes 


Yes 


Yes 


0.8 


6.1 


0.2231 ± 0.0060 


0.1195 



We have got the following values for parameters in parameterizations of parton distributions 
(at Ql = 20 GeV2): 





= 1.44, 




= 2.06, 


^ns 


= 1.87 




= 3.88, 




= 3.84, 


'^ns 


= 4.23, 


"'NS 


= 10.9, 


"ns 


= 4.04, 


"'NS 


= 5.03, 



(48) 



where the symbols P, D and F denote the parameters of parameterizations for proton, deuteron 
and iron data, respectively. 

We note that the values of the coefficients are close to ones obtained in other numerical 
analyses (see |jT2|, ^ and references therein). The values of fejy^g (/ = P,D,F) are in 

quite good agreement with quark-counting rules of Ref.||72|]. There is also good agreement with 
theoretical studies ||73|, ^ . 



Table 5. The values of the twist-four terms. 





hJr ) nf Hn 

IL4yXi J UL 122 


/t4l^Xj J Ui 1^2 




ZL oLdL 




(J. 25 


-0.149 ± 0.015 


-0.1/6 ± 0.014 


u.oo 


-u.ioi m u.uio 


-u.i(o m u.uiz 


0.45 


-0.214 ± 0.012 


-0.147 ± 0.022 


0.55 


-0.228 ± 0.022 


-0.065 ± 0.037 


0.65 


0.024 ± 0.070 


0.053 ± 0.080 


0.75 


0.227 ± 0.154 


0.130 ± 0.131 



The values of parameters of twist-four term are given in the Table 5. We would like to note 
that the twist-four terms for H2 and D2 data coincide with each other with errors taken into 
account. It is in full agreement with analogous analysis [0]. 

We obtain the following results (at x^(F2) = 250, x^{slope) = 6.1 on 7 points): 

a, (20 GeV^) = 0.2231 ± 0.0060 (stat) ± 0.0075 (syst) + 0.0030 (norm) 

(49) 

a,(M|) = 0.1195 ± 0.0017 (stat) ± 0.0022 (syst) + O.OOlO(norm) 

The last error (±0.0010 to a<j(M|)) comes from fits with free and fixed normalizations between 
different data of SLAG, NM and BFP Gollaborations. 

So, the fits of SLAG, NMG and BFP data based on the nonsinglet evolution give for coupling 
constant: 

a^(M|) = 0.1195 ±0.0030 (total experimental error), (50) 
which corresponds to the following value of QGD mass parameter: 

^Ms ~ (^^'i^ ± 37 (total experimental error) j MeV, 

A^ = (321 ± 44 (total experimental error) ) MeV, (51) 

where the error connected with the type of normalization of data are included already to sys- 
tematic error. 

Looking at the results obtained in two previous subsections we see good agreement (within 
existing errors) between the values of the coupling constant a<j(M|) obtained separately in the 
fits of BGDMS data and ones in the fits of combine SLAG, NMG and BFP data (see Eqs. (^)- 
(^) and (^)-(0)). Thus, we have possibility to fit together all the data that is the subject of 
the following subsection. 

4.3 SLAC, BCDMS, NMC and BFP data 

We use the following common x-cut: x > 0.25 and Ycut with = 6 (see the Table 1) for the 
BGDMS data. After these cuts have been incorporated, the full set of data is 797 points. The 
starting point of QGD evolution is Ql = 90 GeV^. 



4.3.1. The results of fits. 



Wc verify liere tlie range of applicability of perturbative QCD. To do it, we analyze firstly 
the data without a contribution of twist-four terms, i.e. when F2 = Ff*^'"'^. We do several fits 
using the cut > Ql^^ and increase the value Ql^^ step by step. We observe good agreement 
of the fits with the data when Ql^^ > 10 GeV^ (see the Table 6). 

Later we add the twist-four corrections and fit the data with the usual cut > 1 GeV^. 
We have find very good agreement with the data. Moreover the predictions for as{M^) in both 
above procedures are very similar (see the Table 6 and Fig. 2). 



Table 6. The values of q;s(M|) and different regimes of fits. 



N of 




N of 


HTC 


X'(F2)/D0F 


as (90 GeV^) ± stat 


«s(M|) 


fits 


cut 


points 










1 


1.0 


797 


No 


2.87 


0.1679 ± 0.0007 


0.1128 


2 


2.0 


772 


No 


1.82 


0.1733 ± 0.0007 


0.1151 


3 


3.0 


745 


No 


1.38 


0.1789 ± 0.0009 


0.1175 


4 


4.0 


723 


No 


1.23 


0.1802 ± 0.0009 


0.1180 


5 


5.0 


703 


No 


1.19 


0.1813 ± 0.0011 


0.1185 


6 


6.0 


677 


No 


1.13 


0.1803 ± 0.0013 


0.1189 


7 


7.0 


650 


No 


1.09 


0.1799 ± 0.0016 


0.1179 


8 


8.0 


632 


No 


1.06 


0.1803 ± 0.0019 


0.1181 


9 


9.0 


613 


No 


1.01 


0.1797 ± 0.0023 


0.1178 


10 


10.0 


602 


No 


0.98 


0.1776 ± 0.0022 


0.1170 


11 


11.0 


688 


No 


0.97 


0.1770 ± 0.0024 


0.1167 


12 


12.0 


574 


No 


0.97 


0.1768 ± 0.0025 


0.1167 


13 


1.0 


797 


Yes 


0.97 


0.1785 ± 0.0025 


0.1174 



We have got the following values for parameters in parameterizations of parton distributions 
(at = 90 GeV^): 





= 2.40, 


^NS 


= 2.46, 




= 2.46, 


A^ 


= 1.65, 






= 3.98, 


iD 


= 3.94, 


"ns 


= 4.08, 




= 4.72, 






= 4.85, 




= 2.38, 




= 1.55, 


d^ 
^NS 


= 7.97, 


(52) 



The values are in good agreement with ones presented in previous subsection. Then all 
discussions given there can be applied here. 

Table 7. The values of the twist-four terms. 





hA{xi) of H2 


hA{xi) of D2 




h4{xi) of C and Fe 




± stat 


± stat 




± stat 


0.275 


-0.221 ± 0.010 


-0.226 ± 0.010 


0.250 


-0.118 ± 0.187 


0.350 


-0.252 ± 0.010 


-0.214 ± 0.010 


0.350 


-0.415 ± 0.233 


0.450 


-0.232 ± 0.019 


-0.159 ± 0.020 


0.450 


-0.656 ± 0.494 


0.550 


-0.122 ± 0.360 


-0.058 ± 0.300 






0.650 


-0.159 d= 0.031 


-0.057 ± 0.031 






0.750 


0.040 ± 0.050 


0.020 ± 0.049 
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Figure 2: The values of a<j(M|) and at different Q^-values of data cutes in the fits based 
on nonsinglet evolution. The black points show the analyses of data without twist-four contri- 
butions. The white point corresponds to the case where twist-four contributions were added. 
Only statistical errors are shown. 

The Table 7 contains the value of parameters of the twist-four term. As it was in the pre- 
vious subsection, the twist-four terms for H2 and D2 data coincide with each other with errors 



taken into account that is in agreement with []I4 



So, the analysis of combine SLAG, NMC, BCDMS and BFP data are given the following 
results: 

• When HT corrections are not included and the cut of is 10 GeV^ at the free normal- 
ization 

xVDOF = 0.98 and a,(90 GeV^) = 0.1776 ± 0.0022 (stat). 



a,(M|) = 0.1170 ±0.0009 (stat) 

When HT corrections are included and the cut of is 1 GeV^ 

xVDOF = 0.97 and a,(90 GeV^) = 0.1785 ± 0.0025 (stat). 



(53) 



a,(M|) = 0.1174 ±0.0010 (stat) (54) 



Thus, as it follows from nonsinglet fits of experimental data, perturbative QCD works rather 
well at g2 > 10 GeV2. 

4.3.2. The study of threshold effects. 

Here we would like to study threshold effects in Q^-evolution of SF F2. Note that at NLO 
level in nonsinglet case the coefficient function of F2 and anomalous dimension do not depend 
on the number Uf oi active quarks. Then, the our study of the threshold effects in Q^-evolution 
of SF F2 is exactly equal to the investigation of a role of threshold effects in the QCD coupling 
constant asiQ"^)- 

To study the threshold effects we consider two types of possible thresholds of heavy quarks: 
= 4m J and = rn^f- First type of thresholds has appeared when a heavy quark with 
the mass nif takes a possibility to be born. The second one lies close to the position of 
"Euclidean-reflected" threshold of heavy quarks. It should play a significant role (see |5^) in 
the as((5^)-evolution. 

A. Let thresholds appear at = 4mj. Then we split the range of the data to three 
separate ones: 

• The values are between 1 GeV^ and 10 GeV^, where the number oi active quarks 
is 3. 

• The values are between 10 GeV^ and 80 GeV^, where the number of active quarks 
is 4. 

• The values are above 80 GeV^, where the number of active quarks is 5. 
Table 8. The values of as(M|) and at different regimes of fits. 



of 

fit 


range 




Ql 


of 

points 


X\F2) 


MS 

± stat 
(MeV) 


MS 

± stat 
(MeV) 


Ai^ 

MS 

± stat 
(MeV) 


a.(M|) 
± stat 


1 


1-10 


3 


5 


195 


124 


400 ± 30 


308 ± 26 


220 ± 23 


0.1187 ± 0.0020 


2 


10-80 


4 


20 


455 


471 




291 ± 17 


208 ± 13 


0.1177 ± 0.0012 


3 


80-300 


5 


90 


190 


143 






199 ± 54 


0.1169 ± 0.0040 



The results are shown in Table 8. The average as(M|) value can be calculated and it has 
the following value: 



«s(M|) 



0.1178 ±0.0010 (stat) 



(55) 



B. Let thresholds appear at Q'l = rni. Then we split the range of the data to two separate 



ones: 



• The values are between 2.5 GeV^ and 20.5 GeV^, where the number Uj of active 
quarks is 4. 

• The Q"^ values are above 20.5 GeV^, where the number n/ of active quarks is 5. 



Table 9. The values of q;s(M|) and ^ at different regimes of fits. 



N of 






Ql 


of 




MS 


MS 


«.(M|) 


fit 


range 






points 




± Stat 
(MeV) 


± Stat 
(MeV) 


± stat 


1 


2.5-20.5 


4 


10 


241 


197 


298 ± 10 


213 ± 8 


0.1181 ± 0.0007 


2 


20.5-300 


5 


90 


558 


533 




187 ± 16 


0.1159 ± 0.0014 



The results are shown in Table 9. The average q;s(M|) value can be calculated and it has 
the following value: 

q;,(M|) = 0.1176 ±0.0006 (stat) (56) 

Thus, we do not find a strong dependence on exact value of thresholds of heavy quarks. 
The theoretical uncertainties due to threshold effects can be estimated for as{M\) ^ 0.0002. 

4.4 The results of the analyses based on nonsinglet evolution 

Thus, using the analyses based on NS evolution of the SLAG, NMC, BCDMS and BFP exper- 
imental data for SF F2 we obtain for as{M\) the following expressions: 

1. When we switch off the twist-four corrections, and put the cut > 10 GeV^, we have 
got at x^/DOF = 0.98 

q;,(M|) = 0.1170 ± 0.0009 (stat) ± 0.0019 (syst) ± 0.0010 (norm) (57) 

or 

as{Ml) = 0.1170 ±0.0023 (total experimental error) (58) 

2. When we add the twist-four corrections, and put the cut > 1 GeV^, we have got at 
xVDOF = 0.97 

q;,(M|) = 0.1174 ± 0.0007 (stat) ± 0.0021 (syst) ± 0.0005 (norm) (59) 

or 

q;,(M|) = 0.1174 ±0.0022 (total experimental error) (60) 

Looking at the results obtained in the Section we see that the central vahic of the coupling 
constant q;s(M|) obtained in the fits (based on NS evolution) of combine SLAG, BCDMS, NM 
and BFP data lies between the central values of the coupling constants obtained separately 
in the fits of BGDMS data and in ones of SLAG, BGDMS, NM and BFP data. All obtained 
values of q;s(M|) are in good agreement within existing statistical errors. 




Figure 3: The study of systimatics at different i/cut values in the fits based on combine siglet 
and nonsinglet evolution. The QCD analysis of BCDMS C^^, H2, D2 data (the case of combine 
evolution): no a Xcut, Qq = 20 GeV^. Thresholds of c and b quarks are chosen at = 9 GeV^ 
and = 80 GeV^, respectively. The inner (outer) error-bars show statistical (systematic) 
errors. 



5 Results of fits of F2: the combined nonsinglet and 
singlet evolution 

At this case, the quite low x experimental data lie at low range and we choose Qq = 20 
GeV2. We use also Nmax — 8- 

The study of the A'max-dependence of the results in the combine nonsinglet and singlet 
case of evolution has been found in Note here only that the analysis in |jl9| shows the 



A^max-independence of the obtained results starting already with N„ 



7. 



5.1 BCDMS ^ff^^i)^ data 

As in the previous Section, we start our analyses with the experimental data []^, H, Q obtained 
by BCDMS muon scattering experiment. The full set of data is 762 points. The starting point 
of QCD evolution is Ql = 20 GeV^. 

As in the nonsinglet evolution case we have studied influence of the experimental systematic 
errors on the results of the QCD analysis as a function of Ycuts, Ycuu and Ycutb apphed to the 
data. Here we use also several sets A^ of the values for the cuts at 0.5 < x < 0.8, which are 
given in the Table 10. 

The absolute differences between the values of for the distorted and undistorted sets of 
data are given in Table 11 and the Fig. 3 as the total systematic error of estimated in 
quadratures. The number of the experimental points and the value of for the undistorted 
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Figure 4: The values of the slope d{\n F2) / d{\n Q'^) at Q"^ = 20 GeV^. The white points 
correspond to the theoretical predictions based on combine singlet and nonsinglet evolution. 
The black points show BCDMS C^^, H2 and D2 data without a Ycut- 



set of F2 are also given in the Table 11 and the Fig. 3. 
Table 10. The values of Ycuts, YcutA and Ycutb- 



N 





1 


2 


3 


4 


5 


Ycut'?, 





0.14 


0.16 


0.18 


0.22 


0.23 


YcutA 





0.16 


0.18 


0.20 


0.23 


0.24 


Ycutb 





0.20 


0.20 


0.22 


0.24 


0.25 



Table 11. The values of q;s(M|) at different values of N. 



N 


number 


X\F2)/D0F 


a, (20 GeV^) ± stat 


full 




of points 






syst. error 





762 


1.22 


0.1992 ± 0.0034 


0.0122 


1 


649 


1.06 


0.2116 ± 0.0042 


0.0096 


2 


640 


1.07 


0.2126 ± 0.0044 


0.0088 


3 


627 


1.05 


0.2152 ± 0.0045 


0.0080 


4 


596 


1.04 


0.2172 ± 0.0047 


0.0076 


5 


590 


1.04 


0.2160 ± 0.0047 


0.0068 
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Figure 5: Notation as in Fig. 4 with one exception: the Ycut with = 5 is taken into account. 



From the Table 11 and the Fig. 3 we can see that the values are obtained for = 1 5 
of Ycut3, YcutA and Fcuis are very stable and statistically consistent. The case iV = 5 reduces the 
systematic error in by factor 1.8 and increases the value of as-, while increasing the statistical 
error on the 27%. 

The importance of the Y-cvX can be shown also in the Figs. 4 and 5, where the slope 
d{\nF2) / dilnQ"^) has been shown at = 20 GeV^. As we can see, there is an essential in- 
provement (the corresponding (slope) decreases in half), when the F-cut has been taken into 
account. 

After the cuts have been implemented (in this Section below we use the set A^ = 5), we have 
590 points in the analysis. Fitting them in agreement with the same procedure considered in 
the Section 3, we obtain the following results: 

a,(20 GeV^) = 0.2160 ± 0.0047 (stat) ± 0.0068 (syst) ± 0.0031 (norm) 

(61) 

a,(M|) = 0.1175 ± 0.0014 (stat) ± 0.0020 (syst) ± 0.0011 (norm) 

As in the nonsinglet case the last error (±0.0011 to as(M|)) comes from difference in fits with 
free and fixed normalizations of BCDMS data |^, |^, |^ having different values of energy. 

So, for the fits of BCDMS data [0, H, || based on complete singlet and nonsinglet evolution 
with minimization of systematic errors, we have the following results (total experimental error 
is squared root of sum of squares of statistical error, systematic one and error of normalization): 



as(M|) = 0.1175 ± 0.0026 (total experimental error) 



(62) 
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Figure 6: The values of the slope d{ln F2) / d{ln Q"^) at = 20 GeV^. The white points 
correspond to the theoretical predictions based on the twist-two approximation of perturbative 
QCD and combine singlet and nonsinglet evolution. The black points show SLAG, NMC and 
BFP experimental data. 



The value of as(M|) corresponds to the following value of QCD mass parameter: 

= (290 ± 20 (stat) ± 29 (syst)) MeV, 

= (2O6 ± 17 (stat) ± 24 (syst)) MeV, (63) 

where the errors connected with the type of normalization of data are included already to 
systematic error. 



5.2 SLAC and NMC H2 + D2 data and BFP Fe data 

We continue our analyses with experimental data [|, ^, |TU| obtained by SLAC, NM and BFP 
Collaborations. The full set of data is 719 points (with the cut > 1 GeV^): 364 ones of 
SLAC, 300 ones of NMC and 55 ones of BFP. The starting point of QCD evolution is Ql = 20 
GeV2. 

As in previous Section we give an illustration of importance of 1/Q^ corrections. First of 
all, we analyze the data applying only perturbative QCD part of SF F2, i.e. Fj"'^. Later, we 
add the 1/Q^ corrections: firstly, target mass ones and later twist-four ones. As it is possible 
to see in the Table 12 and Figs. 6-8, we have the very bad fit (x^ (slope) /DOF = 40), when 
we work only with twist-two part ^2"'^. The agreement with the data is better essentially 
{x^{slope)/ DO F ^ 9) when target mass corrections have been added. The incorporation of 
twist-four corrections leads to very good fit of the data: (slope) /DOF ^ 1.05 (see the Table 
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Figure 7: Notation as in Fig. 6 with one exception: the target mass corrections are taken into 
account for theoretical predictions. 



12 and the Fig. 8) . We note that the statistical and systematic errors are combined in 
quadratures. 

Thus, wc sec that (slope) decreases in 38 times when the corrections has been 
taken into account. 

Table 12. The values of as(M|) and at different regimes of fits. 



N of 


TMC 


HTC 


syst. 


X\F^)/D0F 


{slope) 


a, (20 GeV^) 


«s(M|) 


fits 






error 




for 23 points 


± stat 




1 


No 


No 


Yes 


5.5 


800 


0.2400 ± 0.0017 


0.1241 


2 


Yes 


No 


Yes 


2.2 


179 


0.2153 ± 0.0018 


0.1174 


3 


Yes 


Yes 


Yes 


0.85 


21 


0.2138 ± 0.0058 


0.1170 



Looking at the results in the Table 12, we see the following results for coupling constants 
a,(20 GeV^) = 0.2138 ± 0.0058 (stat) + ±0.0075 (syst) + 0.0030(norm) 

(64) 

a,{Ml) = 0.1170 ± 0.0016 (stat) + ±0.0021 (syst) + O.OOll(norm) 

As in the nonsinglet evolution fits, the last error ±0.0011 to as{M'^) comes from fits with free 
and fixed normalizations between different data of SLAG, NM and BFP GoUaborations. 
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Figure 8: Notation as in Fig. 6 with one exception: the target mass and twist-four corrections 
are taken into account for theoretical predictions. 



We would like to compare the results in the Table 12 with the results of the analyses of 
the data when an additional W'^-cut is taken into account. The inclusion of the W^^-cut is 



very popular (see and references therein) and we fit considering data with variation of the 
W^-cut (and with the standard cut > 1 GeV^). The results of the fits (without twist-four 
correction) are presented in the Table 13 (the systematic errors of the data are included in the 
fits). 



Table 13. The values of q;s(M|) and in fits with different values of ly^-cut. 



of 
fits 


cut 


X\F2)/D0F 


a, (20 GeV^) 
± stat 


(MeV) 
± stat 


Ag^ (MeV) 
± stat 


«.(M|) 
± stat 


1 


2.0 


1.30 


0.2407 ± 0.0013 


400 ± 6 


296 ± 4 


0.1243 ± 0.0004 


2 


4.0 


1.00 


0.2135 ± 0.0018 


280 ± 7 


194 ± 5 


0.1169 ± 0.0004 


3 


6.0 


1.00 


0.2070 ± 0.0023 


253 ± 9 


178 ± 7 


0.1150 ± 0.0007 


4 


8.0 


0.91 


0.2128 ± 0.0043 


277 ± 18 


197 ± 14 


0.1167 ± 0.0012 


5 


10 


0.91 


0.2107 ± 0.0053 


268 ± 22 


190 ± 18 


0.1162 ± 0.0015 



As we can see from the Table 12 (the last line, where twist-four corrections are incorporated) 
and Table 13, the results for as(M|) are in very good agreement for values of W^-cnt larger 
then 4 GeV^. So, the W^-cnt procedure can be used successfully to switch off the range of 
experimental data where higher-twist corrections are required. 



We would like to note that the results obtained in two previous subsections show very good 
agreement between the values of coupling constant a<j(M|) obtained separately in the fits of 
BCDMS data and ones in the fits of combine SLAC, NMC and BFP data (see Eqs. (llD-dlD 
and (|64D). Thus, as in the case of nonsinglet evolution we have possibility to fit togather all 
the data. It is the subject of the following subsection. 

5.3 SLAC, BCDMS, NMC and BFP data 

Here we start to analyze the maximal number of experimental points which have been produced 
in considered experiments. The full set of data is 1309 points. 

5.3.1. The study of range, where l/Q"^ corrections are important. 

Here we would like to repeat our analysis given in Subsection 4.3 . Firstly we fit the data 
without a contribution of twist-four terms. We use the cut > Qlut ^'^d increase the value 
Qcut step by step. Later we do fits including the twist-four corrections and the cut > 1 
GeV^. 

As it was in nonsinglet case, we observe very good agreement for first type of the fits starting 
with Qlut — 15 GeV^ (see the Table 14 and the Fig. 9). For the second type of fits the agreement 
is good already at Q"^ > 1 GeV^. The both types of the fits give very similar results. Moreover, 
the results are very close to ones obtained earlier in the nonsinglet case (see the Tables 3 and 6). 



Table 14. The values of a;s(M|) and at different regimes of fits. 



of 




A^ of 


HTC 


X'(F2)/D0F 


a,(20 GeV^) 


AllL 

MS 


«s(M|) 


fits 


cut 


points 






± stat 


(MeV) 


± stat 


1 


1.0 


1309 


No 


1.55 


0.2258 ± 0.0011 


333 


0.1203 ± 0.0004 


2 


4.0 


1051 


No 


1.27 


0.2364 ± 0.0017 


380 


0.1232 ± 0.0005 


3 


6.0 


942 


No 


1.30 


0.2385 ± 0.0022 


390 


0.1237 ± 0.0005 


4 


8.0 


870 


No 


1.32 


0.2232 ± 0.0035 


321 


0.1196 ± 0.0010 


5 


10.0 


817 


No 


1.27 


0.2226 ± 0.0035 


318 


0.1194 ± 0.0011 


6 


11.0 


793 


No 


1.21 


0.2187 ± 0.0038 


301 


0.1183 ± 0.0011 


7 


12.0 


758 


No 


1.18 


0.2192 ± 0.0039 


304 


0.1185 ± 0.0011 


8 


13.0 


754 


No 


1.17 


0.2180 ± 0.0039 


297 


0.1181 ± 0.0012 


9 


14.0 


740 


No 


1.17 


0.2169 ± 0.0041 


294 


0.1178 ± 0.0013 


10 


15.0 


714 


No 


1.14 


0.2177 ± 0.0042 


297 


0.1180 ± 0.0013 


11 


1.0 


1309 


Yes 


1.11 


0.2167 ± 0.0024 


293 


0.1177 ± 0.0007 



We obtain the following results: 

• When twist-four corrections are not included and the cut of is 15 GeV^ 

xVDOF = 1.14 and a,(20 GeV^) = 0.2177 ± 0.0042 (stat), 

a^(M|) = 0.1180 ± 0.0013 (stat) (65) 

• When twist-four corrections are included and the cut of is 1 GeV^ 

xVDOF = 1.11 and a,(20 GeV^) = 0.2167 ± 0.0024 (stat), 

a,(M|) = 0.1177 ±0.0007 (stat) (66) 
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Figure 9: The values of ^^(Ml) and at different Q^-values of data cutes in the fits based on 
combine singlet and nonsinglet evolution. The black points show the analyses of data without 
twist-four contributions. The white point corresponds to the case where twist-four contributions 
were added. Only statistical errors are shown. 



For additional illustration of importance of l/Q^ corrections at nonlarge values we study 
the slope d{ln F2) / d{lnQ'^) as it has been done in the previous subsection 5.2 . First of all, we 
analyze the data applying only perturbative QCD approximation of SF F2 (with target mass 
corrections taken into account), i.e. Ff*^*"^. Later, we add the cut > 15 GcV^. As it is 
possible to see in the Figs. 10 and 11, we have the bad fit {x^ (slope) / DO F ^ 7.78) in the case 
without a Q'^ cut. The agreement with the data is strongly better when this cut has been 
added: (slope) ^ 1.26 in the case. 

As in the previous subsection the incorporation of twist-four corrections leads also to very 
good fit of the data (without a cut): (slope) /DOF ^ 1.09 (see the Fig. 12) . These 
results demonstrate the importance of twist-four corrections at nonlarge values. 

Thus, as it follows from the fits of experimental data based on combine singlet and nonsinglet 
evolution, perturbative QCD works well at > 15 GeV^. 



5.3.2. The study of threshold effects. 
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Figure 10: Tlie values of tlie slope d{ln F2) / d{ln Q'^) at = 20 GeV^. Tlie white points 
correspond to the theoretical predictions based on perturbative QCD (with target mass correc- 
tions taken into account) and combine singlet and nonsinglet evolution. The black points show 
SLAG, BCDMS, NMC and BFP experimental data without a Q"^ cut. 



Here we continue our study of threshold effects in Q^-evolution of SF F2. Note that at LO 
level and NLO one in the singlet case of evolution the coefficient functions of F2 and anomalous 
dimensions depend on the number n/ of active quarks. 

By analogy with the NS case of evolution (see subsection 4.3.2), to study the threshold 
effects we consider two types of possible thresholds of heavy quarks: = 4m j and Qj- = rnj. 
First type of thresholds has appeared when a heavy quark with the mass rrif takes a possibility 
to be born (in the framework of photon-gluon fusion process, for example). The second one 
lies close to the position of "Euclidean-reflected" threshold of heavy quarks. It should play a 



significant role (see [^) in the ^^(Q )-evolution. 



A. Let thresholds appear at = Am'j. Then we split the range of the data to three 
separate ones (see page 20). 



Table 15. The values of as(M|) and at different regimes of fits. 



of 

fit 


range 




Ql 


of 

points 




MS 

± stat 
(MeV) 


A^ 

MS 

± stat 
(MeV) 


A^ 

MS 

± stat 
(MeV) 


a,(M|) 
± stat 


1 


1-10 


3 


3.0 


467 


290 


331 ± 24 


250 ± 20 


176 ± 16 


0.1148 ± 0.0015 


2 


10-80 


4 


20 


627 


595 




274 ± 21 


194 ± 17 


0.1165 ± 0.0014 


3 


80-300 


5 


90 


190 


156 






220 ± 70 


0.1187 ± 0.0050 
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Figure 11: Notation as in Fig. 10 with one exception: the cut > 15 GeV^ is taken into 
account for experimental data. 



The results are shown in Table 15. The average q;s(M|) value can be calculated and it has 
the following value: 



a^(M|) = 0.1158 ± 0.0010 (stat) 



(67) 



B. Let thresholds appear at = mj. Then we split the range of the data to two separate 
ones (see page 20). 

Table 16. The values of q;s(M|) and different regimes of fits. 



TV of 
fit 


range 




Ql 


N of 
points 




MS 

± stat 

(MeV) 


MS 

± stat 

(MoV) 


«3(M|) 

± stat 


1 
2 


2.5-20.5 
20.5-300 


4 
5 


10 
90 


519 
631 


396 
670 


230 ± 21 


160 ± 16 
205 ± 15 


0.1132 ± 0.0016 
0.1174 ± 0.0013 



The results are shown in Table 16. The average q;s(M|) value can be calculated and it has 
the following value: 



as{Ml) = 0.1157 ±0.0020 (stat) 



(68) 



The results are very surprising. From one side, all variables: the coefficient functions of F2 
and anomalous dimensions, depend on the number n/ of active quarks. However, we do not find 
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Figure 12: Notation as in Fig. 10 with one exception: the twist-four corrections are taken into 
account for theoretical predictions. 

a strong dependence on exact value of thresholds of heavy quarks. From another side, the cen- 
tral values of the average as(M|) obtained here are essentially lower than in other our analyses. 

Thus, the theoretical uncertainties due to threshold effects can be estimated in the case of 
combine singlet and nonsing let evolution for as(M|) as 0.0001. 

5.3.3. The values of fitted parameters. 

We have got the following values for parameters in parameterizations of parton distributions 
(at Ql = 20 GeV2) ^ 



au(20) 


= 0.72, 


bum 


= 3.72, 




ad(20) 


= 0.69, 


&d(20) 


= 5.81, 




Cs(20) 


= 0.375, 


bsm 


= 13.8, 




Pom 


= 0.519, 


bom 


= 11.4, 




irf(20) 


= 1.222, 




= 0.554, 


^3^(20) = 0.253, 


K[{20) 


= 1.10, 




= -0.081, 


7^3^(20) = -0.58 



(69) 

For the coefficients 0^(20) and ad(20) we find good agreement between their values and the 
double- logarithmic estimations in Refs. ||75|, [76[] , based on We would like to note that the 



estimations in Ref. |7g] have been given in other set of parameters that changes effectively only 



^°Here and in the following subsection we give the results for the coefficient Pg{Qo) but not for the one 
{Qq). They are connected because of Eq. (pi| ) 
beta-function B{a, b) has been defined in Eq.(p6^ 



Cg{QI)- They are connected because of Eq. (^ ): Pg(Qo) = C'g((9o) ' B{ag{Ql) + l,bg{Ql) + 1), where the 




Figure 13: The values of the nuclear-effect ratio: r"^, for A = C^"^ and Fe^^. The solid 
and dashen curvess correspond to the small-x asymptotics ~ x~'^ of sea quark and gluon 
distributions with = and uj = 0.18, respectively. 



the value of normalization point Q^. As it has been shown in Refs. 0, the value of a„(20) 



and arf(20) should be nearly Q^-independent (if the values are not too close to 1) ^ This 
Q^-independence of values of au(20) and 0^(20) explains our good agreement with the results 
of |7^. The values of a„(20) and 0^(20) are supported also by recent fits (see discussions in 



Ref. 



The value &u(20) is in agreement with Eqs.(^) and (0) and with other fits p!2| , pT^ , |26| , 



that supports its slow Q -dependence (see ||73|, ^). The value of 6^(20) is higher than 6„(20) 
that is supported by other fits (see, for example, [0, |13|) and references therein) and by quark- 
counting rules The values of 6g(20) and 65(20) are very high, that is in agreement with 
BCDMS analyses [0, ||, || and demonstrates difficulties to study the large-x asymptotics of sea 
quark and gluon distributions in analyses of inclusive deep-inelastic data 0. 

The value of Pci^O) shows that at Qq = 20 GeV^ gluons contain about half of nucleon 
momentum. 

The coefficients i^f (20) and (20) {i = 1, 2, 3) demonstrate non-zero values of nuclear 
effects for bound nucleons in and Fe nuclei. The Eqs. (|39|) with the values of coefficients 

^^This Q^-independence is very similar to corresponding Q^-independence of the coefficients 0^(20) and 
00(20) in the power-like small x asymptotics ~ x°'^ and ^ a;""^ of singlet parton distributions, if as and ac are 
not close nmnerically to (see studies in Refs.|^, ^ and references therein). 

In semi-inclusive case of deep- inelastic scattering the gluons give large contributions, essentially at low x 
values, (see, for example, the recent study of open charm production in psf and references therein) and, thus, 
gluon distribution can be perfectly extracted. 
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Figure 14: The values of the twist-four terms. The black and white points correspond to the 
small-x asymptotics ~ a;""^ of sea quark and gluon distributions with uj = and uj = 0.18, 
respectively. The statistical errors are displayed only. 



K[^{20) and Kf {20) given in Eqs. ( |69D demonstrate the shapes of nuclear effects which are 
represented in Fig 3, where we see a resonable agreement of our curves with the experimental 
data from Refs fS^, 

The values of twist-four terms are given in the Table 17. To obtain the values we used the 
approximate equality of twist-four terms for H2 and D2 targets that has been obtained in our 
studies in the previous Section (see the Tables 5 and 7). This is also in agreement with the 
Ref. The values of twist-four terms are represented also in Fig. 14. 

We would like to note (see the Table 17 and Fig. 14) about a quite strong rise of twist- 
four terms at lower x-bins. The necessity of large magnitude of twist-four corrections at the 
low X values it is possible to observe also in the Figs. 6, 7 and 10, where there is a quite 
strong difference between experimental data and theoretical predictions (based on perturbative 
QCD) for the slope d{lnF2)/d{lnQ'^). The rise is in good agreement with theoretical predic- 



tions 1^ and with the recent analyses of HI and ZEUS data at low values of x and (see P3|). 



Table 17. The values of the twist-four terms. 





1 o Lex L 






IT ■ 




U.UUo 


0.87 ± O.lo 


r\ r\r\f\ 

0.090 


o.lo ± 0.03 


0.275 


-0.19 ± 0.01 


U.Ulo 


0.8o ± 0.12 


r\ -1 r\r\ 

0.100 


n f\c\ _i_ n no 

0.09 ± 0.02 


0.350 


-0.19 ± 0.01 


U.Uio 


0. /^O ± 0.10 


0.110 


n ntx _i_ n no 
O.OO ± 0.03 


n /I n 


n 1 o _i_ n no 
-U.12 ± U.02 


u.uzo 


u.uo m u.uo 


n 1 AD 




u. ouu 




0.035 


0.57 ± 0.06 


0.150 


0.43 ± 0.11 


0.550 


0.04 ± 0.03 


0.050 


0.39 ± 0.04 


0.180 


-0.13 ± 0.01 


0.650 


0.35 ± 0.05 


0.070 


0.28 ± 0.03 


0.225 


-0.15 ± 0.01 


0.750 


0.66 ± 0.10 


0.080 


0.30 ± 0.15 


0.250 


-0.27 ± 0.13 







5.3.4. BFKL-like parameterizations of gluon and sea quark distributions. 

As we have already discussed in Section 2, we would like to try to study the parameters of 
the sea quark and gluon distributions when the terms ~ x°'^^'^o) and ~ x"''-^^'^^^ were incorpo- 
rated. These terms take into account a possible rise of the sea quark and gluon distributions at 



low X values. As it has been already noted in Section 2, from DGLAP-like analyses |4^, ^ 
the parameters as and ac should be the same, because they are mixed together into the "-|-"- 
component of the Q^-evolution (see [0). Moreover, the parameter u = —as = — og should be 



Q -independent (see, for example, ||80|, ^), if it is not small, i.e. x >> Const at small x. 

In the fit with free nonzero co value we have got the following values for parameters in 
parameterizations of parton distributions (at Ql = 20 GeV^) 0: 



a„(20) 


= 0.72, 


bum 


= 3.69, 




ad(20) 


= 0.68, 




= 5.44, 




a5(20) 


= -0.18, 


acm 


= -0.18 




C5(20) 


= 0.185, 


bsm 


= 10.4, 




Pg(20) 


= 0.524, 


bom 


= 7.31, 




i^f(20) 


= 1.160, 


K^m 


= 0.472, 


K^{20) = 0.141, 


i^f (20) 


= 1.03, 


K^m 


= 0.131, 


K^{20) = -0.2 



(70) 

We would like to note that the values of parameters of valent quark distributions are not 
changed really. The values of &g(20) and 65(20) are yet high but they are closer to predictions 
of quark-counting rules |7^ than the corresponding values obtained in the previous subsection. 



The values of the parameters of the nuclei effect ratio are not changed within considered 
errors. The similarity of the results for the nuclei effect ratio is shown in the Fig. 13. 

The value of uj is equal to 0.18, that is in perfect agreement with the recent studies based on 
BFKL dynamics [^] when NLO corrections [^, ^ were taken into account (see, for example, 
studies a review and references therein). Moreover, this value is in good agreement 
also with recent phenomenological studies (see a recent review in |^0[) of Pomeron intercept 
values and also with recent HI and L3 data ^ . 



We would like to note that the fit contains strong correlations between the values of ui, the coupling 
constant and twist-four terms. These correlations come because of very limited numbers of experimental data 
used here lie at the low x region. Indeed, only the NMC experimental data contribute there. Then, the results 
( [zoj ) can be considered seriously only when HI and ZEUS data |43[ Q have been taken into account. We hope 
to incorporate the HERA data 0, H in our future investigations. 



As it is possible to see in the Tables 17 and 18 and also in Fig. 14, the effect of strong 
rise of twist-four magnitude at small x values observed in previous subsection is completely 
absent here ^ So, the rise is replaced by the small x rise of twist-two gluon and sea quark 
distributions. This replacement seems due to a small number of experimental points at low x 
range and narrow range of values there. The cancellation of twist-four corrections at low x 
is in good agreement with the recent studies 92|. This demonstrates the fact that a strong 



rise of twist-four corrections coming from BFKL-like approaches [|8^ has negligible magnitude 



Table 18. The values of the twist-four terms. 





hi{xi) 


Xi 

± stat 


hi{xi) 
± stat 


Xi 


hi{xi) 
± stat 


0.008 


0.004 ± 0.090 


0.090 


0.11 ± 0.03 


0.275 


-0.14 ± 0.02 


0.013 


0.05 ± 0.09 


0.100 


0.05 ± 0.02 


0.350 


-0.17 ± 0.02 


0.018 


0.09 ± 0.09 


0.110 


0.05 ± 0.03 


0.450 


-0.12 ± 0.03 


0.025 


0.11 ± 0.08 


0.140 


-0.01 ± 0.02 


0.500 


0.43 ± 0.23 


0.035 


0.13 ± 0.07 


0.150 


0.62 ± 0.12 


0.550 


0.01 ± 0.05 


0.050 


0.11 ± 0.05 


0.180 


-0.07 ± 0.02 


0.650 


0.26 ± 0.08 


0.070 


0.11 ± 0.04 


0.225 


-0.09 ± 0.02 


0.750 


0.47 ± 0.12 


0.080 


0.31 ± 0.16 


0.250 


-0.16 ± 0.14 







The value of q;^5(M|) in the fit (with the number of points = 1309 and /DOF = 1.1) 
is as follows: 

a,(M|) = 0.1187 ±0.0015 (stat), (71) 

i.e. it is in good agreement within statistical errors with fits performed earlier but the middle 
value is slightly higher. 

5.4 The results of analyses with combine singlet and nonsinglet evo- 
lution 

Thus, using singlet analyses of the SLAC, NMC, BCDMS and BFP experimental data for SF 
F2 we obtain for as(M|) the following expression: 

a,(20 GeV^) = 0.2167 ± 0.0024 (stat) ± 0.0080 (syst) ± 0.0012 (norm) 

(72) 

a^(M|) = 0.1177 ± 0.0007 (stat) ± 0.0021 (syst) ± 0.0005 (norm) 



Looking at the results obtained in the Sections we see very good agreement between the 
value of coupling constant as(M|) obtained in the fits of combine SLAC, BCDMS, NMC and 
BFP data and the values of as(M|) obtained separately in the fits of BCDMS data and in ones 
of SLAC, BCDMS, NMC and BFP data. 

As it was in previous subsection, to obtain the values we used the approximate equahty of twist-four terms 
for H2 and D2 targets that have been obtained in our studies in the previous Section (see the Tables 5 and 7). 
This is also in agreement with the Ref. H. 



6 The dependence on factorization and renormalization 
scales 



In the section we study the dependence of our resuhs on the different choice of the factorization 



scale hf and the renormahzation one /jr. Following the studies |Tj, ^ we choose three following 
values (1/2, 1, 2) for the coefficients kp and kji. 

6.1 Nonsinglet evolution case 

The results are given in the Table 19. We do fits here without higher-twist corrections (no HTC), 
with the number of points 596, at > 10.5 GeV^ and for free normalization of different sets 
of data. The change of the value of coupling constant as(M|) at some kp and k^ values is 
denoted by the difference: 

Aa,(M|) = a,(M|)-a,(M|)|fc^=fc^=i (73) 



Table 19. The values of as(M|) at different values of kp and fc^. The values in brackets 
correspond to the case when the Eg . (|2^) replaces the Eq.(|20|) into the NLO corrections to 
coefficient functions. 



kR 


kp. 




a, (90 GeV^) ± stat 


«.(M|) 


A«,(M|) 


1 


1 


556 


0.1789 ± 0.0023 


0.1175 





1/2 


1 


558 


0.1769 ± 0.0022 


0.1167 


-0.0008 








(0.1745) 


(0.1155) 


(-0.0020) 


1 


1/2 


545 


0.1730 ± 0.0021 


0.1150 


-0.0025 


1 


2 


568 


0.1876 ± 0.0025 


0.1211 


+0.0036 


2 


1 


555 


0.1826 ± 0.0025 


0.1191 


+0.0016 








(0.1858) 


(0.1203) 


(+0.0028) 


1/2 


2 


570 


0.1856 ± 0.0026 


0.1203 


+0.0028 








(0.1817) 


(0.1186) 


(+0.0011) 


2 


1/2 


554 


0.1770 ± 0.0022 


0.1167 


-0.0008 








(0.1784) 


(0.1173) 


(-0.0002) 


1/2 


1/2 


556 


0.1789 ± 0.0023 


0.1175 


-0.0034 








(0.1694) 


(0.1134) 


(-0.0041) 


2 


2 


567 


0.1912 ± 0.0028 


0.1225 


+0.0050 








(0.1965) 


(0.1245) 


(+0.0070) 



We find similar variation of as(M|) with the variations of kp and /cr: as{M'^) increases 
(falls) with increasing (decreasing) of values of kp and/or kji. So, the dependence is quite 
similar to one which has been obtained in ^ ^ by the variation of /c-scale from 1/4 to 4 
{k = kp = kR in H, MM)- 

Taking maximal and minimal values (that corresponds to kR = kp = 1/2 and 2, respec- 
tively) of coupling constant we obtain the theoretical uncertainties +0.0050 and —0.0034 for 
as(M|). In the case when the replacement ( p2[ ) has been used also in NLO corrections to the 
coefficient functions (i.e. when the Eq.(p2|) replaces the Eq.([20|) there), the theoretical uncer- 
tainties for as(M|) are little higher: +0.0070 and -0.0041. 



Thus, using the analyses with NS evolution of the SLAC, NMC, BCDMS and BFP experi- 
mental data for SF F2 we obtain for a;s(Mf ) the following expressions (when no HTC, > 10 
GeV^ and = 0-98): 

a,(M|) = 0.1170 ± 0.0009 (stat) ±0.0019 (syst) ±0.0010 (norm) + { QQ34 (theor),(74) 
or 

a,(M|) = 0.1170 ±0.0023 (total experimental error) ± { QQ34 (theor), (75) 

where the symbol theor marks the theoretical uncertainties which contain the sum of the scale 
uncertainties, threshold error (±0.0002) and the method error (±0.0002) in quadratures. 



6.2 Combine singlet and nonsinglet evolution 

The results are given in the Table 20. We do fits with higher- twist corrections, with the number 
of points 1309, at > 1 GeV^ and for free normalization of different sets of data. 



Table 20. The values of as(M|) at different values of kp and kn. The values in brackets 
correspond to the case when the Eg . (|2^) replaces the Eq.(20) into the NLO corrections to 
coefficient functions. 



kR 


kp- 




a, (20 GeV^) ± stat 


MS 

MeV 


Ai^ 

MS 

MeV 


«s(M|) 


Aa,(M|) 


1 


1 


1410 


0.2167 ± 0.0024 


293 


209 


0.1178 





1/2 


1 


1410 


0.2112 ± 0.0019 


270 


191 


0.1162 


-0.0016 






(1443) 


(0.2104 ± 0.0029) 


(267) 


(189) 


(0.1160) 


(-0.0018) 


1 


1/2 


1423 


0.2040 ± 0.0020 


241 


168 


0.1140 


-0.0038 


1 


2 


1447 


0.2300 ± 0.0031 


351 


256 


0.1215 


±0.0037 


2 


1 


1413 


0.2204 ± 0.0024 


309 


222 


0.1189 


±0.0011 






(1500) 


(0.2263 ± 0.0030) 


(334) 


(242) 


(0.1204) 


(±0.0026) 


1/2 


2 


1422 


0.2190 ± 0.0029 


303 


217 


0.1185 


±0.0007 






(1500) 


(0.2132 ± 0.0031) 


(278) 


(197) 


(0.1167) 


(-0.0011) 


2 


1/2 


1460 


0.2021 ± 0.0022 


233 


162 


0.1134 


-0.0044 






(1496) 


(0.2323 ± 0.0030) 


(361) 


(264) 


(0.1220) 


(±0.0042) 


1/2 


1/2 


1436 


0.1975 ± 0.0012 


216 


149 


0.1120 


-0.0058 






(1450) 


(0.1970 ± 0.0018) 


(214) 


(148) 


(0.1120) 


(-0.0058) 


2 


2 


1447 


0.2340 ± 0.0033 


369 


271 


0.1225 


±0.0047 






(1460) 


(0.2343 ± 0.0032) 


(370) 


(271) 


(0.1226) 


(±0.0048) 



We find that variations of q;s(M|) with the variations of kp and kji are very similar to ones 
which have been obtained in previous subsection. However, there is a quite big difference in 
the cases k^ = 2, kp = 1/2 and k^ = 1/2, kp = 2 between results in the Table 20 in brackets 
and without ones. The difference seems to come from the correlations between the values of 
higher-order contributions (that is mimicked by scale dependences) and twist-four corrections, 
i.e. so-called duality effect (see and references therein). 

As in the case of nonsinglet evolution, the dependence of q;s(M|) with the variations of kp 
and ku is quite similar to one which have been obtained in |^ by the variation of /c/j-scale 
from 1/4 to 4. 



Taking maximal and minimal values (that corresponds to = kp = 1/2 and 2, respec- 
tively) of coupling constant we obtain the theoretical errors +0.0050 and —0.0057 for as{M'^). 
In the case when the replacement (^) has been used also in NLO corrections to the coefficient 
functions (i.e. when the Eq.(^) replaces the Eq.(^) there) the theoretical uncertainties for 
as{Mz) are changed very little but x^(-^2) is higher. 

Thus, using these analyses of the SLAG, NMC, BCDMS and BFP experimental data for 
SF F2 we obtain for 

«,(M|) = 0.1177 ± 0.0007 (stat) ± 0.0021 (syst) ± 0.0005 (norm) + | QQ57 (theor), (76) 

where the theoretical uncertainties contain the scale ones (see above), the ones due to threshold 
effects (±0.0001) and the method error (±0.0002) in quadratures. 

In conclusion of the Section we would like to note that the theoretical uncertainties in both 
types of analyses (based on nonsinglet evolution and on combined singlet and nonsinglet one) 
are essentially larger than the corresponding total experimental errors. 



Indeed, the total experimental errors are as follows: 
in the analyses with the nonsinglet evolution: 

(total experimental error) = 

r (stat) + (syst) + (norm) = 0.0038 (total linear experimental error) 

2 2 2 (77) 

(stat) + (syst) + (norm) = 0.0023 (total quadratic experimental error) 

in the analyses with the combined singlet and nonsinglet evolution: 



(total experimental error) 



0.0033 (total linear experimental error) 
0.0023 (total quadratic experimental error) 



i.e. they are less by factor 1.5 2 to compare with corresponding theoretical uncertainties. 

As it has been shown in [^, ^ |2^, the theoretical uncertainties decrease essentially 
(by a factor around 2.5), when NNLO corrections have been taken into account. So, the fits of 
combined data show real necessity in analyses of DIS data at NNLO approximation. 



7 Summary 



As a conclusion, we would like to stress again, that using the Jacobi polynomial expansion 
method, developed in |l^, [1^, we have studied the Q^-evolution of DIS structure function 
F2 fitting all modern experimental data existing at values of Bjorken variable x: x > 10~^. 



1. From the fits we have obtained the value of the normalization as{M'^) of QCD coupling 
constant. First of all, we have reanalyzed the BCDMS data cutting the range with large 
systematic errors. As it is possible to see in subsections 4.1 and 5.1 (and also the Figs. 1 
and 3), the values of as(M|) rise strongly when the cuts of systematics were incorporated. In 
another side, the values of as(M|) does not dependent on the concrete type of the cut within 
modern statistical errors. 

The values as(M|) obtained in various fits are in good agreement with one other. Indeed, 
we have very similar results for as(M|) in separate analyses of BCDMS data (with the cuts of 
systematics) and other ones. This gives us the possibility to fit all data together. 



We have found that at > 10 -r- 15 GeV^ the formulae of pure perturbative QCD (i.e. 
twist-two approximation together with target mass corrections) are in good agreement with all 
data. The results for as{M'^) are very similar for the both types of analyses: ones, based on 
nonsinglet evolution, and ones, based on combined singlet and nonsinglet evolution. They have 
the following form: 

from fits, based on nonsinglet evolution: 

a,(M|) = 0.1170 ± 0.0009 (stat) ± 0.0019 (syst) ± 0.0010 (norm), (79) 

from fits, based on combined singlet and nonsinglet evolution: 

a,(M|) = 0.1180 ±0.0013 (stat) ±0.0021 (syst) ±0.0009 (norm), (80) 

When we have added twist-four corrections, we have very good agreement between QCD 
(i.e. first two coefficients of Wilson expansion) and data starting already with = 1 GeV^, 
where the Wilson expansion should begin to be applicable. The results for as(M|) coincide for 
the both types of analyses: ones, based on nonsinglet evolution, and ones, based on combined 
singlet and nonsinglet evolution. They have the following form: 

from fits, based on nonsinglet evolution: 

a,(M|) = 0.1174 ± 0.0007 (stat) ± 0.0019 (syst) ± 0.0010 (norm), (81) 

from fits, based on combined singlet and nonsinglet evolution: 

a^(M|) = 0.1177 ± 0.0007 (stat) ± 0.0021 (syst) ± 0.0009 (norm), (82) 



Thus, there is very good agreement (see Eqs. (^), (pOD, ( |8TD and (0)) between results 
based on pure perturbative QCD at quite large values (i.e. at > 10 -i- 15 GeV^) and the 
results based on first two twist terms of Wilson expansion (at >1 GeV^ where the Wilson 
expansion should be applicable). 

We would like to note that we have good agreement also with the analysis of combined 
HI and BCDMS data, which has been given by HI Collaboration very recently. The shapes of 
twist-four corrections are very similar to ones from [l^, Our results for as(M|) are in good 
agreement also with the average value for coupling constant, presented in the recent studies 
(see [|^, EO, |65|, |73, p3, P5[ and references therein) and in famous Bethke review . 



2. As the second item of our summary we would like to note about the real importance 
of NNLO corrections in analyses of DIS experimental data. The incorporation of the NNLO 
corrections have been started already several years ago in various ways (see Introduction for 
discussions) . 

The results are based on the studies of the effect of high order corrections, which can be 
estimated from the dependence of our results on factorization scale fip and renormalization 
one fiR. As it has been point out already in the previous Section the value of the theoretical 
uncertainties 0, coming from this dependence of the results for as(M|) (given by Eqs.(|7^ and 
([TBI) for two types of Q^-evolution), are equal to 

^^As it has been already shown the scale choices /i_F = /i_r = 2Q^ and = IJ-r = Q'^/S give the maximal and 
minimal values of as(Mf ) (at the various choices of values kp = 1/2, kp — 2, k^ — 1/2 and kji = 2 separately) 
and, thus, give the basical part of theoretical error. The additional theoretical uncertainties due to our method 
error and choice of threshold points are negligible. 



A«.(M|)|th, 



eo 



+0.0047 
-0.0057 



^3) 



Thus, the theoretical uncertainties are higher essentially than the total experimental error 
^). Similar values of the theoretical error can be found in recent analyses of DIS process (see 



9|, ^) and of e"'"e~-process in As it has been studied recently by van Neerven and 

Vogt the value of theoretical error decreases strongly (by a factor around 2.5) when the 

NNLO corrections have been taken into account. Thus, our fits of combined data performed 
here and also other analyses [^, 95 1 show real necessity to include the NNLO corrections to 
the study of DIS experimental data. 

As it has been noted in Introduction, using partial information about NNLO QCD correc- 
tions several fits of experimental data have been performed 
references therein). 



see 



and 



In order to do the analyses of experimental data in full range of x values, 
it is necessary to know exactly all NNLO QCD corrections. At present three-loop corrections 
to anomalous dimensions of Wilson operators are still unknown. These calculations, which are 
known only for several finite number of fixed Mellin moments pHI, will be performed 



m 



nearest future by using modern approaches (see |100| , |9 
diagrams. 



to evaluate complicated Feynman 



3. At the end of our paper we would like to discuss the contributions of higher twist 
corrections. 

In our study here we have reproduced well-known x-shape of the twist-four corrections at 
the large and intermediate values of Bjorken variable x (see, for example, the Tables 5, 7 and 
17 and also, for example, the results of very popular article [14|). 

We would like to note about a small-x rise of the magnitude of twist-four corrections, when 
we use flat parton distributions at a; ^ 0. The rise is in full agreement with the theoretical 
predictions [^. As we have discussed already in the Section 5, there is a strong correlation 
between the small- a; behavior of twist-four corrections and the type of the corresponding asymp- 
totics of the leading-twist parton distributions. The possibility to have a singular type of the 
asymptotics leads (in our fits) to the appearance of the rise of sea quark and gluon distributions 
at low X values. At this case the rise of the magnitude of twist-four corrections is 



as 



X 



-0.18 



completely canceled. This cancellation is in full agreement with theoretical and phenomeno- 
logical studies and low x experimental data of HI and L3 Collaborations (see discussions in 
subsection 5.3.4). 



We would like also to give a few words concerning the IRR-model predictions for the twist- 
four and twist-six corrections. 

In our previous study |21|] based on the IRR-model predictions for higher twist corrections, 
we have found a strong correlations between these corrections and the value of coupling con- 
stant. The a;s(M|) value tends to be very small: as(M|) = 0.103 ± 0.002 (stat). This study 
has been supported by fits of DELPHI Collaboration (see ||101|| ) and by some other analyses 

. There is, however, a disagreement with the results of the paper ||102[ , where the twist-four 



corrections in the framework of the IRR-model do not lead to decrease the as(M|) value. In 
our opinion, the situation is not so clear here and it needs more investigations. We hope to 
return to this problem in our future studies. 
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